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There is still significant debate about the symmetry of the order parameter in the heavy-fermion supercon- 
ductor UPd2Al3, with proposals for cos (£3), cos (2A/3), sin(fc3), and e'§ sin^). Here we analyze the tunneling 
spectroscopy of this compound and demonstrate that the experimental results by Jourdan et al 1 are inconsistent 
with the last two order parameters, which are expected to show zero-bias conductance peaks. We propose a 
definitive tunneling experiment to distinguish between the first two order parameters. 

PACS numbers: 



I. INTRODUCTION 

Superconductivity in UPd2Al3 with a transition temperature T c = 
2 K was discovered in 1991 by Geibel et al 2 -, and since that time 
much experimental and theoretical work has been performed, with 
a particular aim of establishing the order parameter symmetry. In 
general, this determination is a crucial first step in understanding the 
pairing mechanism of a given superconductor, so that large efforts 
are generally expended in this direction. 

Evidence for unconventional or nodal superconductivity in 
UPd2Al3 has emerged from a variety of experiments. Feyerherm 
et al 3 measured the Knight shift in UPd2Al3 and found a substantial 
reduction below T c , indicative of singlet pairing. In NMR experi- 
ments Tou 4 and Matsuda 5 found low-temperature nuclear spin lat- 
tice relaxation rate T^ 1 power law (T 3 ) behavior, and the absence 
of a Hebel-Slichter coherence peak below T f . Both behaviors are 
characteristic of nodal superconductivity, and in particular the obser- 
vation of T 3 behavior is strongly suggestive of line nodes in the order 
parameter. Similar evidence was obtained by Sate 6 , who measured 
the low-temperature specific heat and found T 2 behavior, also indica- 
tive of line nodes. Hiroi et al 7 measured the low-temperature thermal 
conductivity and also found T 2 behavior. More recently, Watanabe et 
al** measured the angle-dependent magnetothermal conductivity and 
found evidence for a single line node parallel to the basal plane of 
the hexagonal UPd2Al3. 

Tunnelling spectroscopy can be a powerful probe of order param- 
eter symmetry. The pioneering work in this field was performed by 
Blonder, Tinkham and Klapwjik (BTK) 9 , who provided a simple so- 
lution to the problem of Andreev reflection— in s-wave supercon- 
ductors in an N-I-S contact, and by Tanaka et alii, who was able 
to explain the zero-bias conductance peaks (ZBCP) generally ob- 
served in tunneling measurements of the high-temperature cuprate 
superconductors^. Since that time much theoretical work on tunnel- 
ing spectroscopy in various materials has been performed; primary 
references of interest are the book by Tinkhami^ and the study by 
Honerkamp and Sigristi 4 .. 

The basic conclusion of Tanaka's original work was that zero-bias 
conductance peaks are intimately tied in with the phase difference 
between the pair potentials of the transmitted electron-like and hole- 
like quasiparticles on the superconducting side of the junction. When 
this phase difference reaches 31, Andreev reflection (resulting in the 
transmission of two electrons for a single incoming electron) is en- 
hanced and normal reflection diminished, so that the transmission 
coefficient is increased. This enhancement is most prominent at zero- 
energy. At increasing tunneling energy barrier height (the thin oxide 
layer between the superconductor and normal) the normal state con- 
ductance decreases, so that the relative conductance, or ratio of su- 
perconducting to normal state conductance, increases. It should be 



noted that while strongly enhanced ZBCPs (diverging in the limit of 
large barrier height) only result from a % phase difference, signifi- 
cant zero-energy states can be realized in other cases, as ini 4 .. We 
show below that for the case of interest, substantial amounts of zero- 
energy states can be realized even for order parameters without the 
requisite n phase change, when the effects of quasiparticle lifetime 
are accounted for. The magnitude of these zero-energy states, as well 
as the height and energy of the quasiparticle coherence peaks, are 
strong functions of the order parameter, and for the geometry appro- 
priate to UPd2Al3, of the normal metal used in the junction, so that 
a wide range of interesting and informative experimental results can 
be obtained, offering a definitive resolution to the question of order 
parameter symmetry in UPd2Al3. 



II. MODEL 



In general, the pair state of a superconductor is described by the 
Bogoliubov-deGennes equation s 9 ] 11,15 i 16 : 
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Here, as in Klapwijk et al 16 , we take the f's as electron wave func- 
tions and the g's as hole wave functions, and note that the solutions 
to these equations can be written as 
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where u and v are the well-known BCS coherence factors 14 - 16 ] 17 : 
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We note that in the initial BCS formulation these u and v coeffi- 
cients are only defined for \E\ > |A(k)|, but in this work we will 
extend this usage to energies E < A(k) including negative energy, 
and will model the effects of finite quasiparticle lifetime by letting 
E — » E — iT 18 , where F is the quasiparticle scattering rate. Note that 
this produces complex coherence factors whose squared amplitudes 
do not, in general, sum to unity. This is a consequence of the fact that 
the transmitted waves in the superconductor are evanescent waves 
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for E < A(k) and decay rapidly, with the electrons transmitted to the 
condensate 9 . These waves are effectively surface states. This ex- 
tension of the coherence factors is crucial for studying the effects of 
tunneling voltages V < A(k), as are measured. 

With the coherence factors in place, we are now in a position to 
consider the process whereby an incident electron in a normal metal 
undergoes Andreev and normal reflection at the N-S interface, with 
an electron-like quasiparticle and a back-scattered hole-like quasi- 
particle transmitted. This process is depicted in Figure 1. Note that 
the arrows indicate the particle group velocity; for the transmitted 
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FIG. 2: A schematic diagram of the metallic and UPd2Al3 Fermi 
surfaces, for the longitudinal and transverse geometries. For the 
UPd2Al3 Fermi surface we take the main cylindrical sheet 20 and ne- 
glect its corrugation. 



FIG. 1 : A schematic diagram of all the particles involved in the re- 
flection/transmission process at a N-S boundary; the arrows indicate 
the direction of the group velocity. Note that for the holes the mo- 
menta are opposite to this direction. 

hole and Andreev-reflected hole the momentum is opposite to this. 
We have included a delta-function potential H$(x) to model the in- 
clusion of a thin oxide layer between the superconducting UPd2Al3 
and the normal metal. Translational invariance parallel to the inter- 
face dictates that the momentum in this direction be conserved, while 
the momentum perpendicular to the interface is also conserved in the 
approximation where the barrier energy is much less than the Fermi 
energy of the incoming electron. 

We note that unlike in virtually all other studies of the tunneling 
spectroscopy of unconventional superconductors, the Fermi momen- 
tum of the normal metal plays a key role in selecting the wavevectors 
which will contribute to the conductance, as indicated in Figure 2. 
Only a very small portion of the electrons in a spherical Fermi surface 
approximation will lie on the dominating, approximately cylindrical 
Fermi surface of UPd2Al3, as indicated, and therefore contribute to 
the superconducting or normal state conductance. Electrons far from 
the Fermi surface cylinder of the UPd2Al3 will play no role. 

We will also see that the conductance results depend critically on 
whether the wavevector selected lies near a nodal line; thus different 
metals used on the normal side can be expected to display radically 



different behavior. Indeed, we will see that the measured gap magni- 
tude, the zero-energy density-of-states, and even the existence of the 
usual coherence peaks will depend strongly upon the metal used. Fi- 
nally, we have corrected for the effects of the electron effective mass 
mismatch 19 at the interface. 



III. CALCULATION METHODOLOGY 

In our calculations we follow the work of Blonder et al y , Hon- 
erkamp and Sigrist 14 , and Tanaka 11 . The basic method is straight- 
forward: apply the usual boundary conditions at the interface and 
solve for the ampltitudes of the various components. 

The incoming electron wavefunction can be written as 
( )e I * FJCCOse . with the angle between the normal to the interface 
and k. We have suppressed the part of the wavefunction depending 
on the momentum parallel to the interface, as due to the translational 
invariance in this direction all of the involved momenta are equal. On 
the normal metal side, the Andreev-reflected hole, with amplitude a, 
has wavefunction (Ae ilCFXCOS ®, and the normally-reflected electron, 
with amplitude b, has wavefunction ( )e - * y;toose . On the supercon- 
ducting side, the electron-like quasiparticle, with amplitude c, has 
wavefunction f , , ax )e ,l<FXCOse , and the backscattered hole-like 

quasiparticle has wavefunction {^^{^Iff) 1 e ^e'*' rJ ' cos9 . Here the 
exponential factors represent the phase of the gap at the indicated an- 
gle. We are now able to apply the boundary conditions. Continuity of 
the wavefunction at the interface yields the first two of four equations 
for the four unknowns a,b,c and d. 
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The usual boundary condition 9, 1 1 appropriate for a delta-function 
potential H8(x) is 



Vs ,0!-Vv,0; -= ^ffvKO) 



(7) 



However, due to the large effective mass mismatch at the boundary 
(ms ~ lOOmo) we must generalize this condition. Integrating the 
Bogolubov-deGennes equations across the boundary we find 
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with the effective masses as indicated. Substituting the wavefunc- 
tions outlined above into this formula leads to the final two equations 
for the four unknowns a, b, c, and d. We find the following solutions 
for a and b: 



uv 

a(Q,E) = — 



with 



D={u 2 -e i( ^-«-^v 2 )\Z\ 2 m+u 2 



b(e,E) 



>v 2 )(|Z| 2 



D 



(9) 

(10) 
(11) 



Here m = ms/m^ and Z = H/(ivp cos 8) + 1 /2m — 1/2 where vj? is 
the metallic Fermi velocity, and = rmyp i- e - me g a P phase. The 
energy dependence of a and b is contained in the coherence factors u 
and v. We note two interesting effects: 

• The denominator D increases rapidly with mass mismatch 
m and effective barrier height Z. There are therefore narrow 
transmission resonances when the prefactor of the \Z\ 2 m term 
vanishes. For all the order parameters under consideration, the 
factors we an d %-6 are identical, as are ve and \> % -q. When 
the phase factor is unity, i.e. there is no change of order pa- 
rameter sign under the change in angle from to n — 9, the 
prefactor can only vanish when u=v, i.e. E = A(k). Thus 
for the two cosine-containing order parameters, we expectco- 
herence peaks in the differential conductance at the gap en- 
ergy selected by the wavevector matching on the UPd2Al3 
Fermi surface. This implies that the energy, or voltage, of the 
conductance coherence peaks for these order parameters is a 
function of the normal metal's Fermi wavevector, in contrast 
to most tunneling spectroscopy studies, where no effects of 
wavevector matching are considered, because the Fermi sur- 
face on the superconducting side is also considered spheri- 
cal. It is the peculiar geometry engendered by the longitudinal 
Fermi surface of UPd2 AI3 that creates this unusual coherence 
peak effect. 

• For the two sine-containing order parameters, the phase fac- 
tor is -1, so that for a resonance we must satisfy u 2 = — v 2 . 
Strictly speaking this is impossible, since u 2 + v 2 = 1, but the 
factors u and v diverge as \E\ — > 0, so that both the numera- 
tor of a(E, 9) and the second term of D diverge, while the first 
term of D remains finite. As E — > 0, explicit calculation shows 
that, for r = we attain perfect transmission, and therefore, 
we expect zero-energy resonances to appear as zero-bias con- 
ductance peaks for these two order parameters. 

Following BTK and Tanaka 9,1 1 , the normalized conductance at tem- 
perature T = 1/(3 is given by 
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Note that we are not integrating over the full range of angles 6 as the 
majority of such angles will result in wavevectors lying away from 
the nearly cylindrical 20 Fermi surface of UPd2Al3. In practice what 
we have done in all cases is to use the metallic Fermi momentum to 
determine the angle at which this momentum lies on the UPd2Al3 
Fermi surface and then computed k z o (the component of this mo- 
mentum along the Fermi surface). As a rough accounting for actual 
metallic Fermi surface deviations from spherical we have then inte- 
grated over angles corresponding to k, from 0.95 to 1.05 k z Q. For 
alkali metals such as sodium^S this is an overestimation, while for 
other materials - such as tungsten - this understates the Fermi surface 
deviations from spherical. Nevertheless, the results of our calcula- 
tions are not sensitive to the size of the integration window used, 
provided that this window is small compared with the distance be- 
tween the nodal planes. For all calculations we have taken the ef- 
fective mass mismatch m as 100, while the barrier height H/vp has 
been taken as unity. The results are not sensitive to either of these 
parameters. 



IV. ORDER PARAMETER PROPOSALS 

Before presenting the results of our calculations we give a brief 
summary of the order parameters which have been proposed for 
UPd2Al3. Note that in all proposals £3 is taken as perpendicular to 
the ab-plane, and c is the c-axis lattice constant. 

Won et a! 23 proposed A(k) = cos(2c£3) based upon magnetother- 
mal conductivity data 8 -, while McHale et al 2 ^ performed a strong- 
coupling calculation and found that both cosc&3 and sindy had the 
highest transition temperatures, among possible order parameters. 
It was noted 22 that the chiral order parameter exp(«|>) sincA?; where 
exp(i<|)) = k\ + ik2, is also a possible candidate. All of these order 
parameters contain nodal lines perpendicular to the £3 axis, which 
were found in magnetothermal conductivity 8 measurements on this 
material. 

We note that the last two order parameters are odd under reflec- 
tion around the basal plane; for all order parameters this reflection is 
equivalent to replacing 9 in Figs. 1 and 2 and the above equations by 
71 — 9 (note that this is equivalent to reflection around 6 = 7t/2). This 
is immediately apparent for sincA:3, and the geometric diagram below 
shows that it is true for exp(?'(|>) sinckj. Hence, when the tunneling is 
along the £3 axis, this sign change is expected to produce zero-bias 
conductance peaks (ZBCP) 11 in the tunneling conductance. We will 
see that this is indeed the case. 



V. RESULTS 

Below in Fig. 4 is presented the normalized differential conduc- 
tance dl/dV for the four order parameters described above, for the 
case where the tunneling direction is along the c-axis and the normal 
metal is lead, and T <C r c . These are the conditions of the experimen- 
tal data at T = 0.3K by Jourdan et al 1 . We have here taken the quasi- 
particle lifetime broadening factor T to be 54 /jV, which compares 
well with that extracted from de Haas- van Alphen measurements 20 , 
which measured Dingle temperatures of 0. 1 - 0.28 K, corresponding 
to r = 27 — 16fjV . We find excellent agreement with the data of Jour- 
dan et al 1 throughout the whole range modeled for the order parame- 
ters cos(cfc3) and co%(2ck^), while contrary to experiment, a ZBCP is 
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parameters is relevant for UPCI2AI3, although based upon calcula- 
tions by Nishikawa 24 the s'm(ck?,) gap function may have relevance 
for the isostructural superconductor UNi2Al3. This material has in- 
deed been proposed as a triplet superconductor having an odd gap 
function^ 

For completeness, we have computed dl/dV in the transverse di- 
rection, where the tunneling current is in the basal plane, with the 
same parameters as above. The results are shown in Fig. 5. As 
in the first plot, there is little difference between the two cosine or- 
der parameters, while the sin(ck^) result shows virtually no coher- 
ence peak and no ZBCP, but a high zero-energy density- of-states 
(ZDOS). This is a direct result of the selected angle of incidence for 
lead in this geometry falling very near the nodal line for this order pa- 
rameter; this property will be used to full advantage in proposing an 
experiment to distinguish between the two cosine order parameters. 
We have also performed calculations for the cos(cfc3) gap function 



FIG. 3: A schematic diagram of the wavevector geometry for the chi- 
ral order parameter A(k) = e i<s> s'm(ckj). Here k3 is along the vertical 
axis. Note that the transmitted hole-like and electron-like quasiparti- 
cles have the same value of the phase (]), so that no additional phase 
is introduced by the chiral factor and the results are identical to the 
sin(cA/3) case. 



CZ5 
-4— » 

i 



1.5 



1 



> 

^ 0.5 



A data of Jourdan et al (0.3 K) 
cos(k 3 ) prediction 

— ■ sin(k 3 ),sin(k 3 )e"' > prediction 
cos(2k 3 ) prediction 




.5 -0.75 0.75 1.5 
V(mV) 



FIG. 4: The tunneling data of Jourdan at 0.3 K, with lead as the nor- 
mal metal, is compared with predictions for the four order parameters 
shown. 



predicted for the two sine-containing order parameters. Based on this 
data we thus believe it is rather unlikely that either of these two order 
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FIG. 5: The predictions for the ab-plane tunnelling conductance for 
the three order parameters indicated are shown. 

for the entire range of temperatures tested by Jourdan et al, from 
0.3 K to 1.5 K (in this experiment, T c was 1.6 K.). Figure 6 below 
shows excellent agreement with this order parameter over the en- 
tire temperature range, for nearly all energies We stress that this was 
accomplished with a minimum of fitting parameters; the quasiparti- 
cle lifetime broadening rate T was assumed temperature-independent 
(=54/jV) and the maximum energy gap value A(T) for all tempera- 
tures was near the value appropriate for BCS weak-coupling super- 
conductors of this pairing symmetr y 26 ! 27 . Below we have plotted the 
fitted values of A(T) compared to the BCS prediction. . We note that 
based upon the fits of these data, the ratio 2A(0)/T C is approximately 
3.7, indicating that a weak-coupling treatment may be accurate for 
this calculation. For this order parameter, in the weak-coupling limit 
the ratio 2A(0)/T C is 4.3. We also note, however, the slight bump 
in the data of Jourdan in Fig. 4 around V = ImV, which may be 
indicative of strong-coupling to a magnetic-exciton boso n 21 ] 28 . 
We have not performed detailed comparison with experiment for 
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FIG. 6: The fits to the data of Jourdan using the cos(ck^) order pa- 
rameter are shown. 
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FIG. 7: The fitted values of A(T) are shown. 



the cos(2c£3) order parameter, although it is quite possible that sim- 
ilar agreement with experiment could be obtained for this case. We 
believe that this order parameter may be less likely to apply for 
UPd2Al3 for the reasons presented in2&; namely, the formation of 
a resonance peak in inelastic neutron scattering 29 requires, in the 
magnetic exciton scenario, a sign change over points on the Fermi 
surface separated by the antiferromagnetic wave- vector q = ^£3 and 
this does not apply for this order parameter. Nevertheless, we have 
devised what we believe is a definitive experiment to finally deter- 



VI. PROPOSED EXPERIMENT - ORDER PARAMETER 
DETERMINATION 

The proposed experiment makes use of the unusual Fermi surface 
geometry effects present for tunnelling in the c-axis direction. The 
matching of wavevectors in the metal and in UPd2Al3 places sharp 
constraints on which region of the UPd2Al3 Fermi surface is em- 
ployed by the superconducting electrons. In particular, it is possible 
to choose the normal metal so that for the cos(2ck^) order parame- 
ter, the selected wavevector lies very near the nodal line, while for 
the cos(c%) order parameter, the selected wavevector is much more 
distant from the nearest node. We have chosen Calcium for this pur- 
pose as its Fermi wavevector, based on a spherical Fermi surface 
approximation^ lies very near a nodal line of cos(2cA;3) but is far 
from a nodal line of cos(c&3). The same situation would apply for 
Thallium and Cesium, but these are soft and reactive metals^i that 
are not likely to be useful in forming a tunnel junction. 

In Figure 8 we show the results of the calculation for calcium 
at T=0, with the same quasiparticle lifetime broadening factor as 
before. We observe that for the cos(c&3) case, reasonably well- 
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mine the order parameter in UPd2Al3 after fifteen years of study. 



FIG. 8: Predictions for c-axis NS tunnelling with calcium as the nor- 
mal metal are shown. 

developed coherence peaks, at a magnitude of approximately 250/jV, 
and a comparatively low ZDOS (~ 0.25) are present. However, for 
the cos(2cfc3) case, the coherence peaks are much broader and oc- 
cur at significantly lower bias voltages (V ~ 150^V). In addition, 
the ZDOS is much higher - approximately 0.75. These differences 
are sufficiently robust that we believe that a calcium-based tunnel 
junction experiment, if technically feasible, would be sufficient to 
distinguish between these two order parameters. 

For a final experiment of interest, we have computed the T=0 con- 
ductance for the cos(/c3) case for three metals commonly used in tun- 
neling experiments - beryllium, gold, and aluminum - and show the 
results in Figure 9. As in the previous plot, the coherence peak height 
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FIG. 9: Predictions for c-axis NS tunnelling for cos(c&3) with beryl- 
lium, gold and aluminum as the normal metal are shown. 



and energy vary greatly from one metal to another, with concomitant 
variations in the ZDOS. We believe that such experiments would be a 
good validation of the basic model employed in this paper, and would 
provide useful information regarding the assumption that only very 
few states on the UPd2Al3 Fermi surface contribute to the tunneling 
conductance. 



VII. CONCLUSION 



We have analyzed the tunnelling spectroscopy of the heavy- 
fermion superconductor UPd2Al3, and studied four possible gap 
functions for this material: cos(cfc3), cos(2cA:3),sin(cA:3), and 
exp(i§) sin(c&3). We find that the last two gap functions would be 
expected to show zero-bias conductance peaks in c-axis tunnelling 
experiments as performed by Jourdan et al 1 ; these peaks were not 
observed and therefore these order parameters can be excluded. We 
further find that the cos(cfc3) order parameter gives an excellent fit to 
the experimental data of Jourdan at all temperatures, with a minimum 
of fitting parameters and a A(T) appropriate for the weak-coupling 
BCS theory of this material. Finally, we propose a definitive ex- 
periment to distinguish between the cosfVfc}) and cos(2cfc3) order 
parameters. 
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